
function [dc]=source(t,c_in)

dc = c_in*c_in*t^3;

return
end
